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ABSTRACT 

A  new  computer  model,  SHARC,  has  been  developed  by  the  Air  Force  for 
calculating  high-altitude  atmospheric  IR  radiance  and  transmittance  spectra  with  a 
resolution  of  better  than  1  cm*1.  Comprehensive  coverage  of  the  2  to  40  /im  (250  to 
5,000  cm'1)  wavelength  region  is  provided  for  arbitrary  lines  of  sight  in  the  50-300 
km  altitude  regime.  SHARC  accounts  for  the  deviation  from  local  thermodynamic 
equilibrium  (LTE)  in  vibrational  state  populations  by  explicitly  modeling  the  detailed 
production,  loss,  and  energy  transfer  processes  among  the  important  molecular 
vibrational  states.  The  calculated  vibrational  populations  are  found  to  be  similar  to 
those  obtained  from  other  non-LTE  codes.  The  radiation  transport  algorithm  is 
based  on  a  single-line  equivalent  width  approximation  along  with  a  statistical 
correction  for  line  overlap.  This  approach  is  reasonably  accurate  for  most 
applications  and  is  roughly  two  orders  of  magnitude  foster  than  the  traditional  LBL 
methods  which  explicitly  integrate  over  individual  i*ne  shapes.  In  addition  to 
quiescent  atmospheric  processes,  this  model  calculates  the  auroral  production  and 
excitation  of  COj.  NO,  and  NO*  in  localized  regions  of  the  atmosphere.  Illustrative 
comparisons  of  SHARC  predictions  to  other  models  and  to  data  from  the  CIRRIS, 

SPIRE  and  FWI  field  experiments  are  presented. 

1.  INTRODUCTION 

The  calculation  of  infrared  (IR)  radiance  and  transmittance  spectra  is  important  in 
many  areas  of  atmospheric  science,  including  modeling  the  energy  budget,  analyzing  data 
from  remote  sounding  experiments,  and  understanding  molecular  excitation  and  production 
processes.  The  US  Air  Force  has  developed  a  number  of  computer  codes  that  can  address 
these  applications,  including  LOWTRAN,1  MODTRAN,2  and  FASCODE.3  In  the  upper 
altitude  regime,  more  sophisticated  models  are  required  to  describe  deviations  from  local 
thermodynamic  equilibrium  (LTE)  as  radiative  and  collisional  processes  become 
comparable.  An  initial  non-LTE  (NLTE)  computer  model,  HAIRM,4  resulted  from 

^Phillips  Laboratory,  Geophysics  Directorate/OS,  29  Randolph  Rd.,  Hanscom  AFB,  MA  01731-3010 
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analyses  of  the  1977  SPIRE  rocket  measurements.5'7  In  addition,  high  resolution  radiative 
transport  codes,  NLTE8,9  for  the  quiescent  atmosphere  and  AARC10  for  auroral 
conditions,  were  developed.  However,  a  rapid  and  unified  code  spanning  a  wide  range  of 
altitudes  and  conditions  has  been  lacking. 

This  paper  describes  a  new  NLTE  code,  SHARC,11  which  calculates  upper 
atmospheric  IR  radiation  and  transmittance.  Its  combination  of  speed  and  spectral 
resolution  (0.5  cm'1)  should  make  it  useful  for  many  applications.  The  calculational  model 
includes  all  important  IR-active  species  in  the  2-40  /xm  wavelength  region  for  arbitrary 
line-of-sight  (LOS)  paths  between  50  and  300  km.  It  incorporates  the  bands  of  NO,  CO2, 

03,  H20,  OH,  CO,  and  CH4  found  in  the  quiescent  atmosphere,  including  isotopic  bands 
of  C02,  and  H20.  It  also  accounts  for  auroral  production  and  excitation  of  C02,  NO,  and 
NO+  caused  by  the  flux  of  energetic  solar  electrons. 

In  accounting  for  NLTE  effects,  molecular  vibrational  state  populations  are 
calculated  by  explicitly  solving  the  chemical  equations  for  the  excitation  and  relaxation  of 
each  vibrational  state.  As  in  other  NLTE  atmospheric  models  (e.g.,  Wintersteiner  et  al.,9 
Lopez-Puertas  et  al.12,13),  steady-state  kinetics  are  assumed  for  the  quiescent  atmospheric 
processes  of  collisional  excitation,  de-excitation,  energy  transfer,  radiative  decay,  chemical 
production,  and  illumination  by  the  sun,  earth,  and  atmosphere.  The  additional  production 
and  excitation  mechanisms  resulting  from  electron  deposition  during  auroral  storms  are 
included  in  SHARC  using  a  time-dependent  kinetic  model. 

A  profile  of  vibrational  state  populations  is  input  to  the  LOS  spectral  radiance 
module,  which  performs  line-by-line  (LBL)  radiation  transport  calculations  using  an 
equivalent  width  formalism.  This  approach  uses  an  algebraic  approximation  for  the  total 
absorption  by  a  single  isolated  line,  resulting  in  a  considerable  enhancement  in  speed 
compared  to  the  standard  LBL  grid  method,  by  which  we  mean  explicit  integration  over 
individual  line  shapes.  Corrections  are  made  for  line  overlap.  As  will  be  discussed  below, 
typical  radiance  differences  with  this  approach  are  less  than  10%. 

SHARC  is  available  for  use  by  the  scientific  community  and  may  be  obtained  from 
the  Air  Force  Phillips  Laboratory/GPOS,  29  Randolph  Rd.,  Hanscom  AFB,  MA  01731.  ion' 
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The  schematic  in  Fig.  1  illustrates  SHARC’s  module  structure  and  overall  ^ 
calculational  sequence.  The  input  module  is  a  menu-driven  user  interface.  Atmospheric 
temperature  and  species  density  profiles  are  specified  via  an  external  file.  Profiles  are  'bi.tion  / 
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Figure  1.  SHARC  Module  Structure  and  Calculational  Sequence. 

required  for  the  IR-active  species  (NO,  C02,  H20,  O3,  CO,  OH,  and  CH4),  the  major 
atmospheric  species  to  which  they  are  collisionally  coupled  (N2,  02,  O),  and  atomic 
hydrogen,  whose  reaction  with  03  provides  the  main  source  of  OH(v).  Other  input 
parameters  include  LOS  specifications  and,  if  desired,  the  coordinates  of  a  localized 
auroral  region  through  which  the  LOS  may  pass. 

Altitude  profiles  of  excited  vibrational  state  populations  are  calculated  in  the 
chemical  kinetics  and  radiative  transfer  modules  and  saved  in  a  data  file  for  later  use. 
Under  auroral  conditions,  a  time-dependent  chemical  model  then  calculates  the  additional 
production  of  C02,  NO,  and  NO+  arising  from  interactions  of  auroral  electrons.  To 
generate  the  desired  LOS  spectrum,  the  vibrational  state  populations  are  fed  to  the  spectral 
radiance  module,  which  outputs  radiance  and  transmittance  spectra  and  in-band  intensities. 

Condensed  descriptions  of  the  calculations  are  given  below.  Additional  information 
may  be  found  in  the  SHARC  Users  Manual.11  Code  upgrades  in  progress  or  planned  are 
indicated  by  the  dotted  lines  in  Fig.  1.  These  include:  a  generalized  generator  for 
atmospheric  profiles  of  IR-active  species  based  on  scaling  of  MSIS14  profiles;  extending 
the  spectral  region  into  the  near  IR  and  visible  regions;  and  provisions  for  multiple,  distinct 
atmospheric  regions,  for  use  in  modeling  the  solar  terminator,  tidal  and  gravity  waves,  or 
other  atmospheric  inhomogeneities. 
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2.1  Quiescent  Chemical  Kinetics 

The  quiescent  chemical  kinetics  are  handled  primarily  through  separate  reaction  sets 
for  each  IR-active  molecule;  these  sets  are  solved  individually  in  each  layer  of  the 
atmosphere  to  obtain  the  excited  vibrational  state  population  profiles.  Atmospheric  layers 
are  assumed  to  be  horizontally  uniform  and  typically  2  to  10  km  in  height.  The  reaction 
database  riles  contain  the  list  of  vibrational  states  for  each  molecule  and  the  equations  for 
the  chemical  reactions  and  rate  constants  written  in  symbolic  form.  These  equations 
include  chemical  formation  of  excited  vibrational  states,  collisional  deactivation  and 
excitation  (satisfying  detailed  balance),  spontaneous  emission,  and  the  radiative  excitation 
processes  associated  with  absorption  of  radiation  from  the  sun  and  from  the  atmosphere. 
The  size  of  the  database  is  indicated  in  Table  1. 


Table  1.  Quiescent  Radiating  Species  in  SHARC. 


Molecule 

Isotopes 

States 

Reactions 

Bands 

C02 

3 

84 

664 

113 

h2o 

4 

32 

128 

36 

o3 

1 

30 

176 

45 

NO 

1 

3 

5 

3 

CO 

1 

3 

9 

3 

OH 

1 

10 

32 

24 

ch4 

1 

12 

52 

13 

The  SHARC  CHEMKIN  module,  which  is  based  on  the  Sandia  CHEMKIN 
general-purpose  chemical  kinetics  code,15  reads  the  reaction  riles  and  sets  up  the  time- 
dependent  differential  rate  equations.  The  quiescent  vibrational  state  number  density  [M*] 
in  each  atmospheric  layer  is  obtained  from  the  solution  to  the  steady-state  equations.  The 
species  C02(v),  H20(v),  and  OH(v)  are  indirectly  coupled  to  each  other  via  resonant 
energy  transfer  processes  involving  N2(v=l). 16,17  The  steady-state  equations  for  these 
species  are  linearized  by  equating  the  ground  vibrational  state  number  densities  to  the  total 
number  densities  and  then  solving  for  the  N2(v= 1)  population. 

The  collisional  rate  constants  were  obtained  from  the  recent  literature  and  from 
Taylor’s  review18  of  measurements  prior  to  1974;  details  are  given  elsewhere.19  The  solar 
excitation  rates  are  derived  using  the  LOS  radiation  transport  model  discussed  below.  For 
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solar  zenith  angles  greater  than  90°,  MODTRAN2  is  used  to  account  for  the  attenuation  of 
sunlight  by  the  atmosphere  below  50  km.  The  earthshine  excitation  rates  are  expressed  in 
terms  of  an  effective  blackbody  temperature  corresponding  to  the  altitude  where  the 
vibrational  band  becomes  optically  thick  in  a  nadir  view. 

The  calculation  of  the  radiative  excitation  rate,  ra,  of  molecular  vibrational  states 
due  to  emission  originating  from  atmospheric  layers  in  the  50-300  km  altitude  range  is 
performed  in  a  subroutine  dubbed  NEMESIS  and  follows  the  treatment  of  Kumer  and 
co-workers.16,17  For  a  given  vibrational  band,  the  population  enhancements  of  the 
vibrational  states  are  expressed  via  a  set  of  linear  equations, 

=  oft  XJ  VjtM'jl  ■  (D 

j 

where  is  the  probability  for  absorption  of  a  photon  entering  layer  i,  is  the  probability 
that  a  photon  emitted  from  layer  j  will  be  absorbed  in  layer  i,  and  w;  is  the  branching  ratio 
for  re-emission.  Both  the  P  matrix  elements  and  the  otj  depend  on  the  populations  of  the 
lower  vibrational  states.  Therefore,  the  calculation  proceeds  in  stages  to  solve  for 
successively  higher-energy  states.  The  first  CHEMKIN  run  (with  all  ra  =  0)  defines  the 
ground  state  populations  [M;]  and  the  initial  excited  state  populations  [M*J0  for  transitions 
to  the  ground  state.  The  results  are  used  in  NEMESIS  to  evaluate  (P)  and  a,  and  to 
compute  the  radiative  excitation  rates  ra  for  those  transitions.  CHEMKIN  is  then  rerun 
including  the  ra  to  generate  the  corresponding  populations  for  the  next  set  of  vibrational 
transitions.  The  CHEMKIN/NEMESIS  sequence  is  repeated  until  solutions  for  the  highest 
energy  states  are  obtained. 

The  major  calculational  effort  in  NEMESIS  is  computing  the  P  matrix  elements, 
which  involves  a  multidimensional  integral  over  the  location,  direction,  vibration-rotation 
line,  frequency  location  within  the  line,  and  propagation  distance  of  the  emitted  photons 
within  each  layer.  The  integral  is  evaluated  with  the  aid  of  Monte  Carlo  sampling  using 
trial  "photons".  The  calculation  assumes  semi-infinite  plane-parallel  geometry  and  uses  the 
Voigt  line  shape. 

The  CHEMKIN/NEMESIS  calculations  were  verified  by  comparing  vibrational 
temperatures  for  NLTE  vibrational  states  of  C02,  which  are  both  solar  and  earthshine 
pumped,  with  calculations  by  other  workers.  These  profiles  of  vibrational  temperatures  are 
consistent  with  the  results  of  Wintersteiner  et  al.9  and  Lopez-Puertas  et  al.,20  who  used 
different  computer  algorithms  but  similar  chemical  kinetic  models  and  atmospheric 
profiles.  Typical  results  from  SHARC  are  shown  in  Fig.  2.  The  vibrational  states  are 
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Figure  2.  Calculated  Vibrational  Temperatures  for  C02  States  for  a 
Solar  Zenith  Angle  of  82°. 

labeled  using  the  HTTRAN21  convention.  The  fine  structures  in  the  profiles  are  numerical 
noise  from  the  Monte  Carlo  calculations. 

2.2  Auroral  Kinetics  and  Radiation  Transport 

The  auroral  kinetics  model  calculates  the  enhancements  of  C02,  NO,  and  NO+ 
radiation  resulting  from  energy  deposited  in  the  upper  atmosphere  by  solar  electrons. 
Approximately  600  time-  and  energy-dependent  rate  equations  are  used  to  calculate  the 
secondary  electron  distribution  (in  fourteen  energy  bins)  and  the  subsequent  reactive  and 
energy-transfer  processes.  The  energy  deposition  model  for  the  primary  electrons  is  based 
on  work  by  Grun,  Rees,  and  Strickland22'24  as  implemented  in  the  Air  Force  Phillips 
Laboratory  code  AARC.10  The  chemical  reactions  and  energy  transfer  processes  used  in 
SHARC  are  described  in  Refs.  19  and  25. 

The  integration  of  the  time-dependent  differential  equations  is  accomplished  using 
the  Gear  algorithm  supplied  with  the  Sandia  CHEMKIN  code.13  The  use  of  time- 
dependent  kinetics  results  in  an  improvement  over  the  steady-state  treatment  used,  for 
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example,  in  AARC.10  According  to  SHARC,  the  NO(v=l)  production  efficiency  per  ion 
pair  increases  with  time  during  the  aurora.  This  results  from  collisional  excitation  of  the 
increasing  concentration  of  ground  state  NO,  whereas  in  AARC  the  efficiency  remains 
constant. 

The  emissions  from  NO  and  NO+  are  both  optically  thin  and  prompt;  by  prompt, 
we  mean  that  the  excited  states  radiatively  decay  before  they  have  a  chance  to  relax 
collisionally.  Therefore,  coupling  of  the  radiative  transfer  and  collisional  processes  is  not 
required,  and  the  solutions  to  the  excitation  mechanisms  for  each  layer  suffice.  However, 
C02  j>3  emission  below  100  km  is  both  optically  thick  and  delayed,  so  that  coupling  the 
CO2  time-dependent  chemical  kinetics  to  a  radiative  transfer  scheme  is  required.  With  the 
present  kinetic  scheme,  a  full  treatment  would  require  either  solving  approximately  4000 
coupled  differential  equations  or  repeating  the  radiative  excitation  calculations 
(CHEMKIN/NEMESIS)  at  each  time  increment.  Since  these  are  impractical,  a  modified 
escape  function  approximation,25  which  assumes  that  photons  emitted  from  a  given  layer 
either  escape  the  atmosphere  or  are  re-absorbed  within  the  layer,  has  been  adopted.  To 
insure  the  correct  limit  for  negligible  auroral  excitation,  SHARC  includes  a  correction 
which  forces  a  match  to  the  exact  NEMESIS  solution  for  ambient  conditions.  This 
correction  consists  of  a  time-independent  source  or  sink  term  added  to  the  differential 
equation  for  excitation  of  C02  vy 

2.3  Line-of-Sight  Spectral  Radiance  Model 

In  order  to  determine  the  LOS  spectral  radiance,  the  LOS  properties  are  represented 
in  terms  of  a  sequence  of  homogeneous  segments.  The  spectral  radiance  contribution  of  a 
single  emission  line,  j,  to  a  single  segment,  i,  is  modeled  by 

w>  -  ir  [  f i  w«  •  w<i-us  ]  •  <2> 

where  v  is  the  central  frequency  of  the  calculational  spectral  interval  of  width,  Av,  Ry  is 
the  emission  source  function,  Wy  is  the  cumulative  equivalent  width  for  the  LOS  path  from 
the  observer  through  the  i’th  segment,  and  £  is  a  correction  factor  which  accounts  for  line 
overlap.  The  total  path  spectral  radiance  is  determined  by  direct  summation  of  1^  over  all 
segments  and  lines. 

The  emission  source  term  is  derived  directly  from  the  kinetically  determined 
excited-state  vibration  rotation  populations  and  is  given  by 
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where  Cj  is  the  first  radiation  constant  and  -yy  is  the  ratio  of  the  upper-to-lower  state 
vibration-rotation  populations.  This  expression  for  Ry  reduces  to  the  Planck  blackbody 
function  in  the  limit  in  which  all  the  vibration-rotation  state  populations  are  specified  by 
the  local  gas  temperature,  i.e.,  LTE  conditions. 

In  the  limit  of  no  line  overlap  (£  =  1)  Eq.  (2)  is  exact  provided  that  the  equivalent 
widths  are  exactly  determined,  such  as  by  explicit  frequency  integration  over  the  line  shape 
function.  However,  this  is  computationally  time  consuming.  By  using  the  approximations 
described  below,  a  rapid  evaluation  can  be  achieved  with  only  a  modest  sacrifice  of 
accuracy  (about  ±  10%). 

The  equivalent  widths  Wy  for  single  lines  are  calculated  using  the  Rodgers- 
Williams26  approximations  for  Voigt  line  shapes, 


+  wj 


/WLwDV1 

V  Ww  )  \  ’ 


(4) 


where  the  subscripts  V,  D,  L,  and  W  refer  to  the  Voigt,  Doppler,  Lorentz,  and  weak-line 
limits,  and  aD(cm*1)  is  the  Doppler  linewidth.  WD,  WL,  and  Ww  are  calculated  from  the 
approximations  given  by  Ludwig  et  al.,27  which  are  accurate  to  within  8%  or  better,  using 
Curtis-Godson28,29  path-averaged  line  parameters.  The  line  strengths  and  air-broadened 
half-widths  are  tabulated  in  a  file  generated  from  the  HITRAN  atlas21  that  has  been 
supplemented  with  lines  for  NO+  and  higher  vibrational  states  of  NO  and  03. 

The  line  overlap  correction  is  based  on  the  statistical  overlap  approximation  in 
which  it  is  assumed  that  there  is  no  positional  correlation  among  the  lines  in  the 
calculadonal  spectral  interval.  This  approximation  accounts  for  line  wing  overlap  within 
the  spectral  interval  but  does  not  include  the  contribution  of  line  wings  from  lines  centered 
outside  the  interval.  The  correction  factor  is  given  by 

* = Mi  -ri('-i^)]/ Ewu  •  <5> 

j  i 

In  applying  this  spectral  radiance  model,  it  is  implicitly  assumed  that  the 
calculadonal  spectral  interval  width,  Av,  is  larger  than  the  largest  single-line  equivalent 
width  for  the  endre  LOS.  This  limits  the  calculadonal  spectral  resolution  to  about  0.5  cm"1 
for  tangent  paths  near  the  lower  atmosphere  boundary  of  50  km.  However,  higher  spectral 
resolution  can  be  used  for  non-limb  paths  or  limb  paths  with  higher-altitude  tangent  points. 


-8- 


The  accuracy  of  this  LOS  radiance  algorithm  has  been  explored  through 
comparisons  to  standard  high-resolution  LBL  grid  calculations  for  typical  atmospheric 
conditions.  A  number  of  different  regimes  were  investigated,  including  vertical  and 
horizontal  viewing  geometries,  optically  thin  and  thick  lines,  multiple  overlapping  lines 
and  bands,  and  LTE  and  NLTE  conditions.  Figures  3a  and  3b  show  SHARC  limb 
radiance  calculations  for  just  C02  in  the  v2  spectral  region  at  tangent  altitudes  of  SO  and  75 
km,  respectively.  Each  figure  also  includes  the  absolute  value  of  the  difference  spectrum 
which  was  obtained  by  subtracting  a  LBL  grid  method  calculation  from  the  SHARC 
calculations.  A  major  source  of  error  in  the  SHARC  calculations  is  the  statistical  line 
overlap  approximation.  This  approximation  under-corrects  for  line  overlap  when  the 
emission  lines  in  a  given  spectral  interval  are  nearly  degenerate.  This  can  be  seen  in  the  50 
km  limb  calculation  shown  in  Figure  3a  for  the  Q  branch  near  597  cm'1  and  leads  to  a  20% 
over  prediction  for  the  radiance  over  a  2  cm'1  spectral  interval.  The  line  overlap 
approximation  over-corrects  for  line  overlap  when  the  emission  lines  in  a  spectral  interval 
are  nearly  equally  spaced.  This  occurs  in  the  calculation  shown  in  Figure  3a  in  the  Q 
branch  near  667  cm'1  and  leads  to  an  error  of  10%.  For  LOS’s  and  bands  which  have  less 
line  overlap,  such  as  the  75  km  limb  calculation  in  Figure  3b,  peak  errors  remain  below 
10%.  The  SHARC  calculation  is  roughly  two  orders  of  magnitude  faster  than  a  LBL 
calculation  using  FASCODE. 

3.  DATA  COMPARISONS 

While  comparisons  with  other  computer  codes  can  verify  the  numerical  algorithms, 
upper  atmospheric  IR  data  are  needed  to  assess  the  reasonableness  of  the  various  model 
assumptions,  including  kinetic  mechanisms,  rate  constants,  and  properties  of  the 
atmosphere.  To  illustrate  typical  levels  of  agreement  between  field  data  and  model 
calculations  in  diverse  applications,  several  comparisons  of  SHARC  predictions  and  data 
from  US  Air  Force-sponsored  field  experiments  are  presented. 

The  CIRRIS-1A  experiment30  conducted  on  the  Space  Shuttle  STS-39  during  April, 
1991,  collected  the  first  extensive  set  of  limb  IR  radiance  spectra  over  a  wide  altitude  range 
using  a  sensitive  cryogenic  Michelson  interferometer  with  a  spectral  resolution  of  1  cm'1. 
A  typical  quiescent  nighttime  spectrum  in  the  780-1250  cm'1  spectral  region  is  shown  in 
Figs.  4a  and  4b;  the  tangent  height  at  the  center  of  the  detector  FOV  is  about  64  km.  For 
clarity  in  plotting,  the  resolution  of  the  original  spectrum  has  been  degraded  to  3  cm'1 
FWHM.  The  major  features  are  C02  hot  bands  near  791  cm'1  and  961  cm'1,  the  v3 
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Figure  3a.  SHARC  (Solid)  and  SHARC  minus  LBL  Difference  (Dotted) 
Spectra  for  C02  Radiance  for  a  SO  km  Limb  View. 


Figure  3b.  SHARC  (Solid)  and  SHARC  minus  LBL  Difference  (Dotted) 
Spectra  for  C02  Radiance  for  a  75  km  Limb  View. 
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Figure  4.  Quiescent  Nighttime  Limb  Spectrum  from  the  CIRRIS-1A 
Experiment30  Near  64  km  Tangent  Height  and  SHARC 
Simulation. 
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(1042  cm'1)  and  vj  (1103  cm*1)  bands  of  03,  the  v4  band  of  CH4  near  1311  cm'1,  and 
resolved  lines  of  the  v2  band  of  H20,  which  dominate  the  spectrum  beyond  around  1370 
cm'1.  The  overlaid  SHARC  calculation,  performed  for  similar  viewing  conditions, 
reproduces  the  observed  spectral  structure  throughout  this  region.  The  up  to  factor-of-two 
differences  in  the  absolute  intensities  of  some  bands  are  ascribed  to  differences  in  species 
concentrations  between  the  actual  atmosphere  and  the  model  atmosphere,  which  is  derived 
from  the  NRL  Trace  Gas  Climatology  database.31 

Data  for  upper  atmospheric  C02  p3  emission  and  its  dependence  on  solar  angle  is 
provided  in  Fig.  S,  which  shows  the  limb  radiance  measured  in  an  earlier  rocket 
experiment,  SPIRE,32  launched  from  Poker  Flat  Research  Range,  Alaska  near  the  dawn 
terminator.  Different  phenomena  dominate  the  excitation  depending  on  the  altitude 
regimes.  During  the  day,  the  primary  isotope  is  solar-excited  above  110  km,  while  at 
lower  altitudes  the  radiance  is  enhanced  by  pumping  of  the  2.7  /an  band  and  by  emissions 


O 


Figure  S.  C02  4.3  fim  Band  Limb  Radiance  Measured  in  the  SPIRE  Rocket 
Experiment32  and  Predicted  by  SHARC  for  Different  Solar  Zenith 
Angles  Near  the  Dawn  Terminator. 
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from  hot  bands  and  minor  isotopes.6,33  As  the  density  increases,  the  widening  line  shape 
makes  the  single-quantum  excitation  of  the  primary  isotope  again  important  with 
contributions  from  multiple  quantum  bands  at  1.6,  2.0,  and  2.7  urn.6  The  SPIRE  data  and 
the  SHARC  calculations  agree  to  within  the  data  scatter,  which  averages  around  a  factor  of 
two. 

An  example  of  SHARC ’s  capability  to  calculate  auroral  radiance  enhancements  is 
shown  in  Fig.  6.  A  nitric  oxide  spectrum  from  the  Field-Widened  Interferometer  rocket 
experiment,34  which  observed  an  IBC  Class  II  aurora,  is  compared  with  a  SHARC  model 
calculation  at  a  resolution  of  1  cm'1.  The  auroral  calculation  uses  energy  deposition 
parameters  described  by  Picard  et  al.35  The  reaction  of  metastable  nitrogen  atom,  N(2D), 
with  molecular  oxygen  has  been  assumed  to  be  the  only  source  of  NO  chemiluminescence 
in  the  SHARC  auroral  model.  The  overall  band  shape,  which  includes  strong  aurorally 
induced  hot  band  contributions,  is  reproduced  well  by  the  calculation,  except  near  the 
edges  of  the  band  where  NLTE  rotational  populations,  not  currently  modeled,  yield 
enhanced  high  rotational  lines.  The  calculated  absolute  radiance  is  a  factor  of  four  lower 


FREQUENCY  (cm*1) 

Figure  6.  Calculated  and  Observed34  NO  Spectrum  for  Zenith  Viewing 
Through  a  Class  II  Aurora  from  90  km  Altitude. 
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than  the  data.  Solomon36  and  Gerard  et  al.37  have  suggested  that  the  reaction  of 
transladonally  hot  N(4S)  with  02  may  be  an  important  contributor  to  NO  formed  in  the 
thermosphere.  Sharma  et  al.38  have  presented  calculations  that  show  that  the  reaction  of 
N(4S)  atoms  with  02  accounts  for  recent  observations  of  highly  rotationally  excited  NO 
vibrational  emissions  in  the  dayglow.39  Finally,  the  reactions  of  N^S)  and  N(2D)  with  02 
may  lead  to  a  quantitative  explanation  of  NO  formed  in  recent  artificial  auroral 
experiments.49  Given  the  uncertainties  associated  with  characterization  of  the  aurora  and 
as  well  as  the  role  of  N(4S)  in  the  formation  of  NO,  the  agreement  between  the  model  and 
data  is  encouraging. 

Numerous  other  data  and  model  comparisons  have  been  and  are  currently  being 
performed  for  all  the  major  atmospheric  IR  emission  bands  during  daytime,  nighttime,  and 
auroral  conditions.  The  results  will  be  reported  in  future  papers. 

4.  CONCLUSION 

A  new  computer  m^el,  SHARC,  has  been  developed  by  the  Air  Force  for  rapid 
LBL  calculation  of  NLTE  upper  atmospheric  IR  radiance  and  transmittance  spectra  with  a 
resolution  of  1.0  cm'1  or  better.  SHARC  treats  the  important  molecular  vibrational  bands 
from  2  to  40  /xm  (250  to  5,000  cm'1)  for  arbitrary  lines  of  sight  in  the  50-300  km  altitude 
range,  accounting  for  the  detailed  production,  loss,  and  energy  transfer  processes  among 
the  important  vibrational  states.  Calculated  vibrational  temperatures  agree  with  results 
from  other  NLTE  codes,9,10,13  and  the  equivalent-width  spectral  algorithm  used  in  LOS 
radiation  transport  calculations  results  in  a  considerable  time  savings  over  grid  LBL 
methods  that  explicitly  integrate  over  the  full  line  shape. 

Comparisons  of  SHARC  radiance  predictions  with  field  measurements,  especially 
those  from  CIRRIS  1A,  are  ongoing.  Comparisons  performed  to  date  indicate  reasonable 
agreement  for  most  emission  bands,  including  the  C02  4.3  pm  feature,  which  poses  a  good 
test  for  different  aspects  of  the  code.  Planned  upgrades  include  a  solar  terminator  module 
and  extension  of  the  spectral  region  down  to  the  near  IR.  Data  simulations  using  SHARC 
have  applications  to  remote  sensing  of  the  upper  atmosphere,  such  as  deriving  density 
profiles  for  species  like  NO,  03,  and  H20,  which  have  significant  atmospheric  variability. 
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